Self-consistent Monte-Carlo simulation of 
the positive column in Argon 
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A high accurate self-consistent Monte-Carlo method including charged and neutral particle mo- 
tion and resonance photon transport is developed and applied to simulation of positive column 
discharge in Argon. The distribution of power loss across the tube is investigated. It is shown that 
negative slope of current-voltage gradient curve can not be explained by only one factor. Diffusion 
of metastable atoms to the tube wall in the case of low current ( 0.2 — 3 mA) responds for the most 
part of negative slope. For the high currents the variation of radiation power loss makes a main 
contribution to the amplitude of negative slope. 

PACS numbers: 52.80.Hc, 34.80 
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I. INTRODUCTION 

The positive column of gas discharge has been experi- 
mentally studied for a long time but still there are some 
lack of the simulation techniques which are able to repro- 
duce discharge properties in closed details. 

Modeling of the positive column (PC) is based mostly 
on fluid approaches, which are inapplicable for sheath re- 
gion as well for low pressures and small currents because 
they do not reflect the physical picture of discharge (due 
to assumption of Maxwellian electron distribution func- 
tion and simplified treatment of photon transfer) and re- 
sults in inadequate simulation. 

As it was well established, electron component of dis- 
charge generally reaches a state far from the thermo- 
dynamic equilibrium. Hence the temporal and spatial 
distributions of anisothermic plasmas can be described 
only on appropriate microphysical basis. There are two 
quite different approaches. One way consists of the for- 
mulation of a kinetic equation and its approximate solu- 
tion. Until recently, most studies of the electron behavior 
in positive columns based on solution of the Boltzmann 
equation pj. The other approach uses the technique of 
the particle-in-cell simulation coupled with Monte-Carlo 
(MC) event selection. The last method has no disadvan- 
tages of the kinetic Boltzmann method, and it is limited 
only the available computing power 0, [j| . 

Electron and ion motions in low temperature plasmas 
may be studied by MC approach in which the path of 
a single particle is traced by computer simulation. The 
particle is considered to move classically during period 
of free flight under influence of a self-consistent electric 
field. This field is determined by electron and ion density 
according to Poison's law. The free flight is interrupted 
with scattering by neutral atoms, ionization, excitation, 
wall loss and other kind of collisions which must be con- 
sidered in the model of discharge. The free flight length 



and outcome of these events are selected randomly from 
known probability distributions. Transport properties of 
the discharge under investigation are found from the time 
average of the behavior of the individual particle. When 
an collision occurs one of the possible collision mecha- 
nisms must be chosen. The probability of each mecha- 
nism is proportional to its frequency. Then the particle 
velocity is modified according to the physical nature of 
the mechanism chosen. The simulation of a number of 
particles may be carried out simultaneously and the time 
dependence of their behavior can be derived from an en- 
semble average. 

Resonance photon transfer is also very important for 
a quantitative simulation of gas discharge tube as well 
the charged particle motion. A new simple Monte-Carlo 
treatment of radiation imprisonment based on Holstein- 
Biberman theory Q is developed and embedded into the 
model. 

We present a high accurate self-consistent Monte-Carlo 
approach to simulation of the positive column of Argon 
discharge in this short report. 



II. 



PARTICLE TRANSPORT PROCESSES 



At present model the first 14 excited Ar levels are 
involved into consideration. The 15th level of Argon 
Ari5(3d a ) is a combined level that represents great num- 
ber of the higher excited atomic levels as one state. We 
assume that Ar may be ionized and excited or de-excited 
by electron-atom collisions according to general formula 
which can be written as: 
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(1) 
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where A is atomic symbol, index i denotes initial state 
(before collision), j corresponds to the final state, and 
a{i,j, u) is a cross section for given pair and elec- 

tron kinetic energy is u . In the case of i < j Eq. rep- 
resents excitations and ionization of atom, elastic scatter- 
ing in the case of i — j , and super-elastic electron-atom 
collisions for i > j . 
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FIG. 1: Solid line correspond to momentum transfer cross 
section for electrons in Argon Q. Note the Ramsauer mini- 
mum below 1 eV. Dash line denotes electron impact ionization 
cross section from ground state of Ar. 



The electron excitation and ionization cross sections 
are taken from Refs. 0,0- After ionization of atom a new 
electron appears. We assume that both the ejected and 
scattered electrons are redistributed isotropically, and 
both electrons randomly share available energy which is 
given by subtracting the ionization energy for a given en- 
ergy level from the kinetic energy of the incident electron. 
At very low pressure of Argon the excitation of inner 
electron shells of atom and multiple ionization becomes 
more probable due to the higher electron temperature. 
The present model contains very simplified calculation 
of double ionization events and may obtain only quali- 
tative description of discharges at electron temperature 
T e > 10 eV. 

The cross section of the super-elastic collisions 
a'(j,i;u') (reverse reaction to Eq.QJ) is calculated ac- 
cording to the principle of detailed balance : 



ugi<j(i,j;u) = u'gja' \j,i;u') 



(2) 



where j > i and v! = u — (Ej — Ei) > is initial 
electron kinetic energy, g is a degeneracy of the energy 
level. Therefore the inverse cross section 



ir • ■ is 9i ( , . E j - E i 
a (j, i; u ) = — 1 + ■ — 
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a(i,j;u' + (E J -E i )). 



(3) 

In the MC model we consider elastic electron- Ar scat- 
tering by using the published momentum transfer cross 
section [8| (see Fig0. On average an elastic recoil en- 
ergy loss of 2m e /M = 2.73 x 10~ 5 of the incident kinetic 
energy are used in all simulations. Here m e is electron 
mass and M is Argon atom mass. After elastic collision 
an electron is redistributed isotropically in the center-of- 
mass (CM) system. 

It is often supposed in positive column models that 
all ions and electrons reaching the tube wall recombine 
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FIG. 2: Momentum transfer cross section and resonant 
charge transfer cross section for Ar + in parent gas. Lines 
are fitted functions, symbols denote experimental data [l6||. 

one another. This assumption is also used in our model. 
Recombination of ions at wall produces the neutral atoms 
in ground state. The program redistributes uniformly 
these atoms within tube domain to maintain balance of 
the atom number density. 

In contrast to elastic electron-atom collisions where en- 
ergy exchange is limited by factor 2m e /M <C 1 , en- 
ergy exchange in electron-electron collision has the same 
order as energy of colliding electrons and moreover a 
huge electron-electron momentum transfer cross section 
(Coulomb cross section). For low axial electric field and 
not too small current, the electron-electron collisions re- 
sult mainly in energy transfer from cold bulk electrons to 
hot tail electrons which are responsible for particle (ion- 
ization) balance in discharge. So the correct treatment 
of electron-electron collisions is one of the key elements 
of the discharge simulation. 

To describe electron-electron collisions by MC simu- 
lation we are using classical Coulomb cross section |l(l : 



a(v) 
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ln(Ajj/6o) 
,1/2 



>V2 



(4) 



where = (fcT e /47re 2 n e ) i/ ' is Debye length and 

bo = e 2 /fj,v 2 — 2e 2 /mv 2 is Coulomb radius, fj, is the 
reduced mass of colliding electrons, v is their relative 
speed, and n e , T e are the local electron density and tem- 
perature. The energy exchange during particle-particle 
(here electron-electron) scattering depends on the veloc- 
ity vectors of a projectile lst-particle and target 2nd- 
particle according to equations: 

v' : = (miVi + m 2 v 2 + m 2 vn)/(mi + m 2 ) (5) 



(miVi + m 2 v 2 — mivn) / (mi + tni) 



(6) 
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where and V 2 are the velocities of particles in the 
laboratory system after collision. The unknown unit vec- 
tor n is assumed to distribute isotropically. 

Calculation of electron-electron collisional frequency 
takes much larger computational power than electron- 
heavy particles collisional frequency because of target 
electron distribution is unknown a priori. Let us consider 
electron collision frequency !/, for some i -projectile 
electron moving with given velocity Vj in a spatial do- 
main V given by a computational mesh: 



i\ o"(K - Vj|)/ V (vj)dvj (7) 



where subscript j is the number of target electrons, 
iV(v) is local electron distribution function of the do- 
main V and n e — N/V is local electron density of the 
number of electrons N within V . For relative velocity 
Vij = |vi- Vj\ : 

^i(vj) = n e J o-(vij)f v (Vj)dvj = n e (v^ cr(v i j)) i 

(8) 

Using velocity data of target particles Eq. JSJ can be writ- 
ten as follows 

^( v *) = U = ^7 ( 9 ) 



V 



N 



V 



N; 



where Nj is number of accounted targets. To reduce 
computational efforts one may use Nj < JV as shown 
in the last part of Eq.©. 

We assume in our model that the local collisional fre- 
quency v(yi) in Eq. @ depends only on electron kinetic 
energy u — mv 2 /2 instead of vector velocity v . Then 
the collisional frequency v(u k ) for electrons having ki- 
netic energy in the k energy bin (u k ,u k + Aitfc) be- 
comes 



Ion transport in the ion-parent-atom case is mostly 
guided by resonant charge transfer because of its cross 
section is large at low collision energies. We use a simple 
fitted function for the charge transfer cross section [ll], 
[13] given by 



v(u k ) 



V 



N k N 



(10) 

here S(ui,u k ) — 1 if itj £ {u k ,u k + Au k ) , and 
8(ui,u k ) = if Ui $ {u kl u k + Au k ) , N k = 
Y]^ S(ui, u k ) is number of projectile electrons in k -bin. 

The electron-ion and ion-ion elastic collision processes 
are also included to the model in the same manner as 
for electron-electron collision. However according to our 
simulation these processes do not play any significant role 
in the PC. 

It should be note that for arbitrary projectile-target 
pair the relative velocity may be a very small value and 
the corresponding Coulomb radius b in Eq.Q becomes 
larger than Debye length and/or spatial mesh size. In 
these cases the electron-electron collision event is as- 
sumed as a null collision because the long range electron 
interactions are already included in solution of Poisson's 
equation. 



a cx = (A- Bine) 2 



(11) 



Here e [cV] is a kinetic energy of a projectile ion in 
CM system, A and B are the fitting parameters: 
A at = 6.65 [A] , B A r = 0.3 [A], see FigEI Due to the 
large cross section the ions and atoms are practically un- 
deflected in the center-of-mass system. Therefore, after 
the resonant charge transition occurs a new ion will ac- 
cept velocity of the neutral atom which are randomly 
generated according normal distribution function for a 
given gas temperature which is assumed to be constant 
across the tube. 

Polarization scattering of ion by a neutral atom can 
be described by momentum transfer cross section ^| 
(without 2^2 term): 



a p i = nal \J (a/al)I H /e, 



(12) 



where ao = ti 2 /me 2 — 0.5292 • 10 -8 cm is Bohr radius, 
Ijj = e 2 /2ao = 13.6 eV is ionization potential of Hydro- 
gen atom, and e — M'v 2 /2 is kinetic energy of relative 
motion of ion and atom (that is in the CM system). The 
polarizability of Argon atom is equal to ctA r / a o — 11-08 . 
The minimal value of a p i is limited by the corresponding 
gas-kinetic (hard sphere) cross section a gk as shown in 
Fig|2] We assume that the ion is redistributed isotropi- 
cally in CM system according to Eas. (|5l6fl after an elastic 
collision with a random neutral atom at a given gas tem- 
perature T g = 20° C . 

Because the correct calculation of ion motion is crucial 
point for the gas discharge modeling, we carry out the 
validity check of these cross sections by comparing the 
simulated ion drift velocity with experimental one. Re- 
duced mobility of Ar + in argon gas is fio ( Ar + , Ar) = 
1.6 cm 2 V _1 s _1 hence the experimental drift veloc- 

ity for gas pressure 0.261 Torr and electric field 2 V/cm 
is 93.18 m/s which is in good agreement with simulated 
Ar + drift velocity 92.3 m/s. 

We take into account excited atom-atom collision ion- 
ization process as follow Ar* + Ar* — > Aro + Ar + + e , 
where Ar* is any one of a number of excited Ar levels. 
Cross sections of these processes found in the literature 
vary from 10~ 20 m 2 to 100 x 10~ 20 m 2 We assume 

the recommended value 24 x 10~ 20 m 2 15]. The new 
ionized electron gains kinetic energy Ei + Ej — Iat > 
7.86 [eV] and the new Ar ion accepts random velocity 
from a neutral Ar atom. 

The simple treatment of diffusion based on equation of 
continuity is included to the model. The flux of atoms is 
estimated on basis of Fick's law given by 



3d 



-DVn(r), D = D (P /P)(T/T ) a , (13) 
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FIG. 3: The transmission probability density at Voigt pa- 
rameter a = 0.0186. Solid line denotes theoretical calculation 
by Eq.lUjJ. 



here Do = 0.073 cm 2 /s is diffusion coefficient of the 
excited Ar atoms in Argon gas at the normal conditions, 
and a = 1.92 0. 

In case of a large gradient of density the diffusion flux 
jd may be larger than kinetic limit of the flux of particles 
with averaged velocity (V a ) given by 



3k = n(V a )/4 



(14) 



To eliminate the overestimation of diffusion processes in 
Ea. (|13fl we assume the total flux as the following formula: 



3 =3k3d/(jk + 3d) 



(15) 



The charged particle is considered to move classically 
during period between collisions (free flight) under influ- 
ence of an self-consistent electric field within discharge 
tube. This field is determined by electron and ion den- 
sity according to Gauss's law: 



<j> EdS = 4ttQ = 4n J (m - n e )dV 



(16) 



Because we suppose particle densities are symmetric 
about tube axis, the electric field E in the cylindrical 
coordinate system has only the radial component E r (r) . 
It provides us the simple formula 



E r {r) = 2Q(r)/r 



(17) 



where Q(r) is the total charge within radius r per unit 
length of discharge tube. 



III. RESONANCE RADIATIVE TRANSPORT 

Due to the high absorption of the resonance photons 
by atoms in the ground state the correct treatment of ra- 
diative transfer in discharge condition is very important 



for simulation model. We apply the Biberman-Holstein 
theory Q of resonance radiative transport for MC simu- 
lation of traveling photons. 

Let us denote p(r) — fcexp(— kr) the probability den- 
sity of a photon captured within [r, r + dr] and without 
absorption along distance [0, r] . Then the probability of 
the absorption anywhere within [0, R] is J p(r)dr = 
1 — exp(— kR) and probability of the traversing a dis- 
tance R is T(R) = 1 — J* 1 p(r)dr — exp(— kR) . Due to 
line broadening (natural, pressure and Doppler broaden- 
ing) the absorption coefficient k depends on a photon 
frequency uj . 

For the pure Doppler broadening case since the velocity 
distribution of the atoms is a Maxwellian with given gas 
temperature, one can obtain the absorption coefficient as 



k(uj) = kg exp[— (lo/uq — 1) c m/2T, 



(18) 



where fco is absorption coefficient in the maximum of 
the line profile, c is speed of light, m is an atomic 
mass, and T g is a gas temperature in energetic units. 

It can be shown that the averaged probability density 
of a photon capture after traveling distance r is 



exp 



-2x 2 - k re x 



dx 



(19) 



here x = — ,,,„ a;o c/(2T g /m) 1 ^ 2 , and fco is given by 



0J 



fco 



Afj 92 



riA 



8tt 3 / 2 5 i T{2T g /m) 1 / 2 



(20) 



Here r is mean radiative lifetime of an excited level , 
A is wavelength, and ua is number density of atoms. 
We assume that the natural argon gas consists of only 
one isotope with atom mass 39.948 a.u. For given set 
of 15 atomic levels we apply the natural mean life times 
as they are listed in NIST Atomic Transition Probability 
Tables @ (see page 10-88). 

In a general case the emission (or absorption) line 
shape is determined not only by the Doppler broadening 
but natural and pressure (Lorentz) broadening as well. 
The combination of these effects may be given by Voigt 
profile h(a,x) [l7| : 



h(a, x) = — 

IT 



cxp{-y 2 )dy 
a 2 + (x — y) 2 



(21) 



and then the emission profile f(x) and the absorption 
coefficient k(x) are written as 



<p(x) = h(a,x)/y^jr, k(x) = koh(a,x). 



(22) 



Here x was defined above and a is the Voigt parameter 
(ratio of the Lorentz HWHM to Doppler width at 1/e 
maximum) defined below 



An 



47tt(2T/to) 1 /2 



(1 + 2n A T(v AA <JA)) , 



(23) 
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where a a is a resonance broadening cross section which 
is a a — 643 A 2 for 106.67 nm Ar line, and a a = 
2339 A 2 for 104.82 nm Ar line in the pure Argon dis- 
charge. In the program we use the alternate formula for 
resonance broadening [17j given by 



Hvaa^a) = 0.4506 



92 Ag 

91 67T 2 



(24) 



The Voigt parameters for Argon discharge at pressure 
0.261 Torr are a = 0.0172 for A = 104.82 nm, and 
a = 0.00468 for A = 106.67 nm. 

The function h{a, x) from Eq. (|21|l has two asymptotic 
limits for a = and a — ► oo . The pure Doppler broad- 
ening is realized as h{a,x) — > cxp(— x 2 ) when a — ► . 
In the case of a — » oo the Voigt profile tends to Lorentz 
profile as h(a, x) — » a/^/rr(a 2 + x 2 )) . 

Similarly to Ea. (|19(l the transmission probability den- 
sity in the most general form is given by 



p(r) = / ip{x)k{x)e~ k[x)r dx 



(25) 



and can be rewritten for the Voigt profile as follow 

p(kor) = —= / h 2 (a, x) exp [(— koh(a, x)r]dx (26) 
V 71 " J 

To generate a random number according to the prob- 
ability density p{k$r) from the previous Equation we 
develop a new simple random generator 



k r = (X- 2 - 1) /2 
with comparison function 

/(kor) = 2/(1 + 2fc r) 3/2 



(27) 



(28) 



At first step the program generates a trial random num- 
ber kor according to Ea. I|27|l by using a random num- 
ber X uniformly distributed within [0,1] . At second 
step the uniformly distributed within [0, f(kor)] ran- 
dom number Y is picked over. If Y > p(kor) then the 
program rejects the trial number kor and return to the 
first step. Otherwise the program accept the trial num- 
ber. The efficiency of this generator is equal 50% because 
of J f(k r)d(k r) = 2 . Figure|21indicates that the prob- 
ability density generated by random generator Eq. is 
in a excellent agreement with theoretical function from 
Eq.l|26p. The transmission probability density p{k a r) is 
tabulated and stored at first start of the program. 



IV. 



PARTICLE BALANCE THEORY 



To maintain the charge balance in a steady PC an 
electron may produce many excited atoms and only one 
electron-ion pair during the electron mean life time. The 
exited atoms may disappear by radiative decay, in elec- 
tron inelastic or super-elastic collisions, in electron im- 
pact and atom-atom ionization processes, and due to dif- 
fusion to the tube wall. In our model the ions can only 



recombine at the tube surface. It is clear that the total 
balance of each sort of particles must be equal zero in the 
steady positive column. 

As example of simplified theoretical model of PC let 
us consider the electron balance equation which can be 
written as sum of the volume rates: 



^ v ie n e + ^2 k K 







(29) 



where first term corresponds to direct ionization from 
i-level atoms by electron impact with frequency Vi e — 
ki e rii in unit volume, second term denotes atom-atom 
ionization in collisions between excited i-level atom and 
j-level atom, and third term is electron loss rate on the 
tube wall. 

In the simplest theory of DC positive column (Schot- 
tky, 1924) the electron-wall loss frequency v ew is as- 
sumed to be independent on electron density n e and 
equal to v ew — D e (2.405/i?) 2 , where D e — p{T e is 
the ambipolar diffusion coefficient and T e weakly de- 
pends upon the electron density. In the simplest model 
of positive column where the ionization from ground 
state is only taken into consideration the ionization fre- 
quency vo e does not depend on electron density and 
ionization rate voejie is linear term with respect to n e . 
The electron balance equation Ea. (|29|l becomes simply 
voe. = v ew , therefore electron balance does not depend 
on electric current which is proportional to n e . Ioniza- 
tion frequency i>oe is a function of the axial electric field 
E z . Hence E z is independent of electric current and 
only determined by electron wall-loss frequency. 

We may extend Schottky theoretical model by includ- 
ing ionization from one sort of the excited resonance 
atoms only as it have been done in [Ig| : 



ko e n n e + k xe n x n e - v e 







k 0x n n e - k xe n x n e - v x0 n x = 



(30) 



(31) 



where v x o — 1 jr x is radiative decay frequency of x- level 
atom to ground state 0, n x is number density of excited 
atoms. By solving Ea. (|31|l for n x and rewrite Ea, 131|) 
to obtain balance equation 



k 0e n n e 1 + 







(32) 



It is clear that the ionization term in Eq. ll-il'l) is superlin- 
ear in n e , except the case k xe n e 3> v x o (high current). 
For this model we may conclude that the ionization from 
resonance states results in negative slope of V-I charac- 
teristic curve 0|- ^ is interesting to note that for high 
current case and/or for metastable atom with v x q = 
the balance equation (|3^|l gives a flat V-I curve. 

Let us consider in contrast to 01 the same model 
Eas. l|30l31|l but now for metastable x- level atoms with 
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v xQ = D X (2A05/R) 2 
the tube wall: 



as frequency of diffusion loss on 



k 0e n n e + k xe n x n e 



(33) 



kQ X n n e - k xe n x n e - u x0 n x = (34) 
and we get the similar solution given by 



k 0e n n e 1 + 



v x q 



Vxon e = 



For this model in the enough low current case ( k xe n e < 
v x o ) the ionization from metastable states causes a neg- 
ative V-I slope. Because the number of metastable atoms 
is much higher than the number of resonance atoms in 
low current region Eq.J3HJ conforms to the discharge 
much better than Eas. H3UI31l) . Hence we may conclude 
the diffusion of metastable atoms to the tube wall is an 
essential requirement for a negative V-I characteristic of 
the low current positive column. 

The atom-atom ionization should be included in any 
realistic model. Let us insert this atom-atom ionization 



term k xx n x 



into our model Eas. (|33l34|l : 



koeWoTle ~\~ k xe Tl x Tl e + k xx TL x V e 



ko x n n e 



v x $n x = 



(36) 



(37) 



By combination of two equations it is easy to get the final 
electron balance formula: 



koenone + k xe wn z e + k xx vj 2j n 2 e - v, 



= 



(38) 



where 



— (2ko e + k 0x )n 

k X eJl e l /f x Q 



Again for the enough low current case ( k xe n e < v x § ) 
Eq. (l38|l shows two superlinear ionization terms, the first 
is atom-atom ionization and the second is electron im- 
pact ionization from metastable atoms. These super- 
linear terms result in negative slope of V-I curve. But 
with increasing of current the atom-atom ionization 
term k xx w 2 n 2 becomes a sublinear term in n e (al- 
most independent of n e ) , but the impact ionization from 
metastable atoms tends to linear term in n e . It may 
cause a positive slope of V-I curve if k xe < k xx w. For 
further increasing of current we have to take in considera- 
tion the electron impact ionization from resonance states 
(see Ea. (|33l34|) ') and transitions between metastable and 
resonance states. It adds to Eq. (|38|l a new superlinear 
ionization term which have to turn V-I curve from the 
growth to the drop of voltage. Therefore the unknown 
whole electron balance equations have to contain the ion- 
ization rates consisting from the sublinear, linear, and su- 
perlinear terms as we show above. The interplay between 
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FIG. 5: The voltage-pressure characteristic curve of Ar PC 
at current / = 50 mA. See notations at Figl^l 



them results in a slope of V-I curve. We may distinguish 
at least three electric current region in accordance with 
slope sign of the V-I curve. Low current region (ioniza- 
tion from metastable atom k xe n e < v x q ) corresponds 
to the negative slope, the middle current region is for 
the positive slope (due to independence of atom-atom 
ionization rate on n e ), and the high current region cor- 
responds again to the negative slope ( due to superlinear 
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ionization from excited atoms k xe > k xx w ). 

We have to note that this theoretical model is not ap- 
plicable for very low current PC where the electron-wall 
loss frequency v ew depends on electron density n e , 
and the electron temperature significantly grows with 
decreasing of current. The general condition to be sat- 
isfied for maintenance of the column can be written as 
I r r e = 1 , where r e is a mean life time of an electron 
( r e = Tj = r e (Ni/N e ) at steady conditions). It is clear 
(and easy to see on simulation data) that lower electron 
density /current discharge results in more strong radial 
electric field E r in bulk plasma and relatively smaller 
one in sheath than in higher current case. Under the 
influence of this radial field the ions are accelerated to- 
wards tube wall much faster in the lower current PC. The 
cause of that the mean ion velocity Vi — /liE — DiVni/rii 
is mostly determined by first term and almost does not 
depends on a variation of ion density profile. Therefore 
the mean life time of ions/electrons is longer for higher 
current PC and, hence, ionization rate per a electron 
must decrease with increasing of current to sustain the 
discharge at balance. This conclusion is quite general 
but it can not directly determine the slope of voltage- 
current curve. Nevertheless, for the simplest model (or 
for enough small current) where the ionization rate de- 
pends strongly on axial electric field E z , the equilibrium 
axial electric field E z have to be also less for higher cur- 
rent positive column |Tgj . 

Finally we may conclude that the developed particle 
balance theory itself can not predict precisely the sign 
of V-I slope because the axial electric field E z responds 
first of all for energy balance in discharge. Energy bal- 
ance equations are much more complicated and detail 
information about energy transfers and conversions can 
be obtained by simulation only. 




radial position, R (mm) 

FIG. 6: The spatial dependency of power loss at pressure 
p = 0.261 Tor and current / = 50 mA. 



TABLE I: Contributions of the partial power dissipations in 
negative V-I slope AE„/ AE Z (%) at gas pressure 0.261 Tor 



Current, mA 


Diffusion 


Ionization 


Radiation 


ion-Ar col. 


0.02 - .05 


13.5 


22.5 


30.0 


27.8 


0.05-0.1 


22.9 


18.6 


31.5 


23.2 


0.1-0.2 


29.2 


16.2 


31.9 


20.4 


0.2-0.5 


43.4 


14.2 


24.0 


17.7 


0.5 - 1.0 


54.2 


14.5 


13.7 


18.4 


1.0-3.0 


50.6 


13.1 


21.7 


15.8 


3.0-10 


31.7 


13.3 


40.4 


20.5 


10-30 


5.0 


7.1 


79.3 


10.9 


30 - 100 


0.8 


3.7 


89.3 


7.0 



V. ENERGY BALANCE IN SIMULATED PC 

Electrons acquire kinetic energy from the axial elec- 
tric field E z , from superelastic collisions, from excited 
atom-atom ionization and lose it in the electron- wall and 
electron-atom collisions. Ions gain kinetic energy mostly 
from the radial electric field E r and lose it in the ion- 
atom and ion- wall collisions. Excited atoms in their turn 
gain/lose energy by the electron hits and radiative tran- 
sitions, and lose energy by diffusion to the tube wall. Our 
MC program provides detailed information on these all 
processes included in the model. According to our simu- 
lation the distribution of power loss across the tube is not 
uniform. Figure[S]shows strong nonuniformity in the spa- 
tial dependencies of elastic electron-atom and ion-atom 
power loss in pure argon discharge simulation. Elastic 
(and charge exchange) ion-atom collisions becomes more 
important than elastic electron-atom power loss near the 
tube wall due to strong acceleration of ions by radial 
electric field, especially in the sheath. 

It is evident that the input electric power must be qual 



to the total wasted power at a steady positive column. 
Energy balance equation can be written as follow: 

W = IE Z = W n (40) 

n 

The program calculates independently the input elec- 
tric power W = IE Z and the dissipated powers W n 
for the each energetic n-channel such as light emission, 
electron-atom and ion-atom elastic collisions, ionization, 
diffusion of excited atoms to the tube wall, electron-wall 
and ion- wall losses. Simulation shows that the power bal- 
ance maintains with a good accuracy. To reveal the in- 
fluence of dissipating n-channel on the axial electric field 
E z we may rewrite Ea. (|4U[) as 

E z =^VF„// = ^S„ (41) 

n n 

where E n is a partial axial field which corresponds to 
n-channel of the total power dissipation. By using the 



8 



partial electric field E n we may easily estimate contri- 
bution of each n-channel to the total energy loss in a 
steady PC and, moreover, determine the real atomic rea- 
sons behind variation of measured axial field E z . 

The MC program is able to predict the experimental 
voltage-pressure and voltage-current curve (see Figures 
I4I5H . The negative slope AE Z in V-I curve on Fig^can 
be decompose on partial components AE n of total elec- 
tric field E z according to equality 1 = J2 n AE n /AE z 
from Ea. H41|) . As it is shown in Table I the most consid- 
erable (not for E z but for derivative dE z /dI ) power 
loss channel is diffusion of metastable atoms to the tube 
wall in the case of low current (0.2 — 3 mA). For the 
high currents the variation of radiation power loss makes 



a main contribution to the amplitude of negative slope. 

The Figures 1 1151 also show radiative efficiency of the 
lamp measured as a ratio of radiated power to the total 
dissipated electric power. The maximum of efficiency co- 
incides with the maximum of the voltage-pressure curve 
at pressure around 0.2 Tor and near 20 mA current. 
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